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One of the most disturbing difficulties in tliinking about tlie cosmological constant is that it is not 
stable under radiative corrections. The feedback mechanism proposed in ref. [1] is a dynamical way to 
protect a zero or small cosmological constant against radiative corrections. Hence, while this by itself 
does not solve the cosmological constant problem, it can help solving the problem. In the present 
paper we investigate stability and gravity in this approach and show that the feedback mechanism 
is both classically and quantum mechanically stable and has self-consistent, stable dynamics. 

[HUTP-03/A044] 



I. INTRODUCTION 



The cosmological constant is one of the most outstanding mysteries in contemporary physics [2,3]. There are actually 
two cosmological constant problems. One is the "old" problem: "why is the cosmological constant vanishingly small 
■<!;;j- . compared to the energy scale of, say, the Planck scale?" The ratio of the critical density of the universe today, than 
T-H ■ which energy density of the cosmological constant must be smaller, to the Planckian energy density is as small as 
gince the cosmological constant is not protected against quantum corrections, there is no manifest reason why 
\l ■ it should be so tiny. We now have the second problem: "why is the dark energy roughly of the same order as matter 
', energy density today?" The observation of high-redshift, Type la supernovae is strongly inconsistent with a universe 
■ with vanishing dark energy, whose typical form is a cosmological constant [4-6]. This implies that there should be 
' a small but non-vanishing dark energy and that it should be comparable to matter energy density today. However, 
\^ i since dark energy and matter scale differently as the universe expands, there seems to be no manifest reason why they 
' should be of the same order today. This problem is often called the coincidence problem. 

: Although there are many approaches to the coincidence problem [7-10], all of them seem to have a common 
weakness: they can be spoiled by other contributions to the cosmological constant. This is because it is a priori 
\ assumed that the old problem can be solved. A tiny correction to the unknown solution to the old problem can 
' produce a relatively huge cosmological constant. For example, in the quintessence scenario [7] it is a priori assumed 
that the asymptotic value of the quintessence potential exactly vanishes. If we added a non- vanishing cosmological 
^ ' constant to the potential then the scenario would not work. Moreover, if there are many proposals to solve the 
coincidence problem and if they can coexist then wc will end up with a too much cosmological constant, provided that 
each contribution is additive. Hence, a good solution to the coincidence problem must include a mechanism which 
almost completely cancels other contributions to the cosmological constant. Unfortunately, no phcnomenologically 
^ , acceptable solution with such a mechanism is known. 

' There are also many approaches to the old problem. Among them are quantum cosmology [11-13], membrane 
creation [14-17], a scalar field with a very small mass [18,19], backreaction models [20,21], eternal inflation [22], extra 
dimensions [23], and so on [24,25]. 

We do not consider the quantum cosmology approach based on the Euclidean quantum gravity partly because it 
is not well understood and may be ill-defined. (For example, see refs. [26,27] for challenging arguments.) There is 
another reason why we do not consider this approach. Quantum cosmology is, in a sense, a theory of the beginning, 
or the initial condition, of the universe. For example, the Hartle-Hawking boundary condition [28] predicts a sharp 
peak at a vanishing cosmological constant, but this value should be considered as the effective cosmological constant 
at the beginning of the universe, including all contributions from many fields and quantum corrections. Hence, if 
this prediction is correct then the universe cannot experience inflation nor any phase transitions after its birth. If 
the effective cosmological constant were zero at the beginning of the universe and if a phase transition occurred 
after that then the cosmological constant now would be negative. If we try to avoid this disagreeable conclusion by 
modifying or changing the proposal then the peak must be moved to a non- vanishing cosmological constant. For 
example, the proposals of Linde [29] and Vilenkin [30] prefer a large cosmological constant. This indicates the onset 
of a Planck scale inflation, but does not seem to provide a solution to the cosmological constant problem without 
anthropic consideration [31]. 

A proposal by Rubakov [24] might relax the cosmological constant in the Brans-Dicke frame to a small value if 
the cosmological constant in the Einstein frame is exactly zero. However, for his proposal to work, it is essential to 
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inchidc an inflaton sector in the Einstein frame. Hence, we anyway need to explain why the cosmological constant in 
the Einstein frame is exactly zero today. A similar proposal by Hebecker and Wetterich [25] not only has a similar 
problem but also predicts a too much deviation from Einstein theory for weak gravity and, thus, is phenomenologically 
unacceptable. 

Dolgov's backreaction model [20] and its variants [21] use the kinetic energy of a running scalar field to cancel the 
bare cosmological constant. In these models the effective gravitational coupling vanishes and, thus, the deviation from 
Einstein gravity is too large. 

As far as the author knows, all other previously known approaches to the old cosmological constant problem 
arc cither dependent of the anthropic principle or phenomenologically unacceptable (eg. a non-standard Fricdmann 
equation at low energy, a non- Einstein weak gravity, a naked singularity, etc.). This situation has led to an increase in 
speculation for the necessity of the anthropic principle [32] . Before resorting to this, it is worthwhile to ask whether 
anything could possibly do what the cosmological constant data requires [3]. It seems likely that the correct way to 
interpret the tiny value of the cosmological constant is that conventional quantum field theory is not the whole story, 
so it is worth seeking acceptable modifications. 

In ref. [1] a new dynamical approach to the cosmological constant problem was proposed. Although this idea by 
itself does not explain why the cosmological constant is so small, it does provide a way to protect a zero or small 
cosmological constant against radiative corrections (without using supersymmetry) and thus can help solving the 
cosmological constant problem. In this approach the true value of the vacuum energy is negative, but the dynamics 
is such that the true minimum is never attained, and the universe settles down to a near zero energy state. 

The purpose of the present paper is to investigate stability and gravity in this dynamical approach. In particular, 
we show that (i) the effective cosmological constant decreases in time and asymptotically approaches zero from above; 
(ii) the evolution of a homogeneous, isotropic universe is described by the standard Friedmann equation at low energy; 
that (iii) classical, linearized gravity in Minkowski background is described by Einstein gravity at distances longer than 
I* = sfoLK and at energies lower than Z^^, where k is the Planck length and a is a dimensionless, positive parameter 
of the model; that (iv) the mechanism is stable under radiative corrections and thus a zero or small cosmological 
constant is protected against radiative corrections; and that (v) quantum fiuctuation at low energy is so small that 
the effective cosmological constant at low energy does not overshoot zero even quantum mechanically. 

This paper is organized as follows. In Sec. II we explain the basic idea in analogy with a simple system of a charged 
particle. In Sec. Ill we review the feedback mechanism proposed in ref. [1] in a more general setting. In Sec. IV we 
derive equations of motion in detail and show that the standard Friedmann equation is recovered at low energy. The 
result in Sec. IV is confirmed numerically in Sec. V. In Sec. VI wc consider reheating in this model. In Sec. VII 
we investigate weak gravity in Minkowski background and show the recovery of the linearized Einstein theory. In 
Sees. VIII we investigate stabilities against radiative corrections and quantum fluctuations. Sec. IX is devoted to a 
summary of this paper. 

II. BASIC IDEA 

One of the most disturbing difficulties in thinking about the cosmological constant is that it is not protected against 
radiative corrections, which usually generate enormous vacuum energies compared to what we observe. The feedback 
mechanism proposed in ref. [1] can be considered as a dynamical way to protect a zero or small cosmological constant 
against radiative corrections. Hence, although the feedback mechanism by itself does not solve the cosmological 
constant problem, it can help solving the problem. In this section we explain the basic idea underlying the feedback 
mechanism. 

Let us consider a simple system of a non-relativistic, charged particle in a background electrostatic potential. The 
action of this simple system is 

(1) 

where V{x) is the electrostatic potential. Let us suppose that eV{x) is an increasing function of x so that the particle 
moves towards the negative x region. Suppose, however, that for some reason, we would like to stop the particle just 
before its crossing x = 0. How can wc stop it? One obvious solution is to control the electrostatic potential so that 
it has a minimum near a; = 0, and introduce a friction. By controlling the potential and the friction carefully, we can 
achieve our purpose. A difficulty in this solution is that there might be forces (eg. strong wind, gravity, someone who 
may kick the particle, etc.) other than the electrostatic force and the friction. These forces make the control of the 
potential extremely subtle and difficult. 




- eV{x) 
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A different, rather radical solution would be available if we could control the mass M of the particle. In this case 
we make M to diverge when x is sufficiently close to x = 0. To make this successful we of course need to watch 
the value of x in real time. This case does not involve a subtle control of anything: x completely stops when M 
diverges, irrespective of the value of eV'{x) and other forces. One might worry that the kinetic term in the action 
might diverge. Actually, this is not the case, and the kinetic term does not diverge but vanishes. Since the equation 
of motion of the position x of the charged particle is 

{Mxy +eV'{x) = Q, (2) 

the momentum tt = Mx remains finite in finite time as far as eV'{x) is finite. (This is the reason why x stops when 
M diverges.) Hence, the kinetic term actually vanishes when M diverges: 

—x^ = — ^ (M ^ oo) 

Indeed, the more singular the mass is, the more quickly the kinetic term vanishes. 

We apply this second idea to a scalar field, whose potential determines the effective cosmological constant but is 
not fine-tuned. The reason why we do not consider a fine-tuned potential, or the analogue of the first idea above, is 
that it is not stable under radiative corrections, which add vacuum energies to the potential and also change the form 
of the potential. We still want to stop the scalar field at or near zero cosmological constant. Let us suppose that a 
scalar field (f) is described by the action 

-^-n<A)], (4) 

where / is chosen so that it vanishes at or near zero cosmological constant. A simple choice might be to make / 
dependent of the value of (/> itself so that / vanishes at a root of the potential. In this case, </> stops when the potential 
vanishes. However, this choice is too naive and would not be stable under radiative corrections: radiative corrections 
add vacuum energies to the potential and, hence, change the potential value at which (j) stops. Indeed, if / is a 
function of (j) only then we can always make a change of variable (j) (j) = J dtj)/-^ so that the kinetic term of the 
new variable ^ is canonically normalized. In this language, the choice of / is nothing but a fine tuning of the potential 
of (f> and, thus, this is never stable against radiative corrections. Hence, we cannot protect zero or small cosmological 
constant against radiative corrections if / depends only on 0. In other words, ^ itself does not know where to stop. 

On the other hand, the spacetime curvature seems to know where to stop (f) since it reflects the value of the effective 
cosmological constant via the Friedmann equation. Hence, instead of considering a (/)-dependent /, we postulate that 
/ depends on the spacetime curvature so that the curvature can tell when should stop. For simplicity, we restrict 
our discussion to the case where / depends on the Ricci scalar R only and vanishes at i? = 0. In this case, cj) stops 
rolling when R = 0. We do not provide a reason why / vanishes at i? = 0. However, we shall show that this choice of 
/ is stable under radiative corrections and leads to a stable dynamics. In other words, instead of providing a solution 
to the cosmological constant problem, wc propose a dynamical mechanism to protect a zero or small cosmological 
constant against radiative corrections (without using supersymmetry) as a first step. The singular kinetic term makes 
the scalar fleld stop at zero curvature, even without a fine-tuned potential. 

With this choice of /, we need to make it sure that R = really corresponds to zero cosmological constant. If the 
kinetic term dominates over the potential term then the value of the potential determined by i? = does not need 
to be zero and, thus, this model does not seem to work. Hence, we are left with the following two possibilities. One 
possibility is that there is a scaling solution to scalar dynamics in which the kinetic and potential energies of the scalar 
field decrease in tandem, so that the field is slowed down when small potential energy is achieved. Unfortunately, in 
our specific realizations, there always seems an instability invalidating the scaling solutions we tried. 

The second possibility, which we shall focus on in the rest of this paper, is that the potential term dominates 
over the kinetic term so that R is determined solely by the contribution from the potential. In this case, the scalar 
field indeed stops at zero cosmological constant automatically. We could make / dependent of (j) as well, but such 
a dependence can be removed from the behavior of the kinetic term near i? = by redefinition of (j) without loss of 
generality. Wc suppose that the kinetic term in (4) represents the term which is the most singular-looking at R = 
among many possible terms in the kinetic part and we did not include less singular-looking terms since they arc not 
important at low energy. We also omitted all other dynamical fields since, as we shall see in the following sections, the 
dynamics of (p is so slow that any dynamical fields other than (f) will settle into their ground state before the universe 
approaches a sufficiently low energy state. The potential V{(f)) includes the ground state energies of all such fields. 

Since / is in the denominator and vanishes in the R limit, the kinetic term threatens to be singular at low 
energy. Moreover, since / is a function of the Ricci scalar R, it introduces corrections to Einstein equation and 
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possibly destabilize the system. On the other hand, from the charged-particle analogue above, we have learned that 
a singular-looking kinetic term is actually regular. The same is true here! The more singular 1// is, the smaller the 
kinetic term and corrections to Einstein equation are! This rather counter-intuitive but very simple fact means that 
by making 1// (for the most singular-looking term) sufficiently singular at i? = 0, we can make corrections to Einstein 
equation as small as we like. In particular, it is indeed possible to make the corrections due to the kinetic term much 
smaller than those from usual higher derivative terms like at low energy so that higher derivative corrections and 
stability are completely controlled by the usual higher derivative terms. Note that if there are several (or many) terms 
in the kinetic part, we do not need to require the singular behavior for all coefficients: just one term is enough. 



III. FEEDBACK MECHANISM 




In the previous section we have learned from a simple analogous system that a scalar field stops in a completely 
regular way when the coefficient of a kinetic term diverges. In this section we consider a scalar field with such a 
singular- looking kinetic term. 

The field Lagrangian we consider is 

R 

2^ + ai?2 + Lkin - V{(l)) 

where / is a function of the Ricci scalar R, k is the Planck length, a and q are constants, and K = —K^d'^(f>d^(f>. Our 
sign convention for the metric is ( — h -|-+). In a word, the idea of this model is to assume a non-standard kinetic 
term. The dependence on curvature is what stalls the field at small cosmological constant. 

We can absorb any nonzero cosmological term into V{(j)) without loss of generality. We assume that (i) the minimum 
of the potential V{(l)) is negative; that (ii) q > 1/2; that (iii) a > 0; and that (iv) the function f{R) is non-negative 
and behaves near i? = as 

f{R) ~ {K^Rf^, (6) 

where m is an integer satisfying 

In the m ^ oo limit, this choice of f{R) may be replaced by f{R) ^ exp[—K~^R~^] or similar functions. We suppose 
that the kinetic term in (5) is the most singular-looking one among many possible terms in the kinetic part. We 
assume that the most singular-looking term behaves like (6) with the condition (7), but other terms do not need to 
behave like that. If there are several terms in the kinetic part then the most singular-looking term is dominant and 
other terms are not important at low energy. Thus, eg. adding the standard canonical kinetic term does not affect 
anything. 

Since / is in the denominator of the kinetic term and vanishes in the k^R limit, the kinetic term threatens to 
be singular at low energy. We shall see below that the large m makes (j) evolve slowly, that the numerator K"^ vanishes 
more quickly than the denominator and that the kinetic term is actually regular. 

Now let us explain the motivations for the assumptions (i)-(iv). (i) We would like to propose a mechanism in which 
the scalar field stops rolling at or near zero vacuum energy. For this purpose, the minimum of the potential Vi^cj)) 
should be negative so that V(0) has a root, (ii) For the stability of inhomogeneous perturbations, it is necessary 
that the sound velocity squared = L^in^K / {'^KLkin,KK + Li.in,K) is positive [33]. In our model this condition is 
reduced to q> 1/2. (iii) The i?-dependence of the kinetic term Lkin produces higher derivative corrections to Einstein 
equation which might destabilize gravity. As we shall explain below, the term aR^ can stabilize gravity at low energy 
if a is positive, (iv) For the term aR? to control the stability, we need to make it sure that the term aR? is dominant 
over Lkin at low energy. As shown below, this is the case if and only if the condition (7) is satisfied. 

We do not know a parent theory that will provide our Lagrangian as the low-energy effective theory. In particular, 
we do not give a reason why / vanishes at i? = 0. In this sense, this model by itself does not solve the cosmological 
constant problem. However, we shall show below that this choice of / is radiatively stable, leads to a stable dynamics 
and could help solve the cosmological constant problem. We treat this model as a purely phenomenological suggestion 
that might motivate further research into the possible parent theory, which is presumably not based entirely on 
conventional four-dimensional field theory. 
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We could make / dependent of (f>, but such a dependence can be removed from the behavior of Lkin near R = 
by redefinition of (j) without loss of generality. Note that Lkin above represents the term which is the most singular- 
looking at i? = among many possible terms in the kinetic part and that we did not include less singular-looking 
terms since they are not important at low energy. We also omitted all other dynamical fields since, as we shall see 
below, the dynamics of (j) is so slow that any dynamical fields other than will settle into their ground state before 
the universe approaches a sufficiently low energy state. 

We now argue that the Lagrangian (5) gives a feedback mechanism that makes the field stall at or near zero vacuum 
energy. The argument for the q = 1 case was alreadv given in ref. [1]. Here, let us simply generalize it to a general q 
(> 1/2). 

In the fiat Friedmann-Robertson- Walker (FRW) background 

ds^ = -dt^ + a{tf{dx^ + dy^ + dz^), (8) 
the equation of motion for a homogeneous (f) is 

7r + 3J?7r + F'(<A)=0, 7r=i-(^^^ , (9) 

where H = d/a, a dot denotes the time derivative and a prime applied to V{4>) denotes the derivative with respect to 
(f). Near V = 0,V can be approximated by a linear function as 

V:^CK-^{cP-<Po), (10) 

where c and (j)o are constants. This approximation is extremely good since, as we shall see, rolls very slowly near 
V = and does not probe the global shape of the potential. 

Without any fine-tuning, the dimensionless constant c should be of order unity. The equation of motion (9) can be 
rewritten as 

where B = k^Htt/c and x = lna{t). This implies that B approaches —(3 — H/H^)~^, if H/H^ changes slowly 
compared to a and is smaller than 3. See Fig. 1 for confirmation of this statement by numerical calculation. Hence, 
the asymptotic behavior of n is 



K^TT ' 



-CK,-^H-^. (12) 



If the kinetic term is small compared to the potential term then the Priedmann equation implies that 

3iJ2 ~ k'^V, (13) 



and (12) can be rewritten as 



Kdt{K^V) ~ CK^^ = cf ■ . (kV)5(5^-2 

~ -(c2)^(kV)'"-5(5?^. (14) 



Hence, we obtain 



(«V)-'^/2~(c2)^LJ^, ^ = 4m-^, (15) 

where to is a constant. From the condition (7), 7 > 3 and 

k'^V^+O {t/K-*oo). (16) 

This result means that the field stalls near V" = 0. See Figs. 4 and 5 for confirmation of this behavior by numerical 
calculation. 
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We have assumed that the kinetic energy is small compared to the potential energy. This assumption is easily 
verified. At low energy {kH <C 1), 

= ~ .-'(.H)'" (-4) ^ 

< k-^{kH)^ < K-^{KHf - y, (17) 

provided that c = 0(1). Here, we have used the condition (7) to obtain the first inequality. We have implicitly 
assumed that the standard Priedmann equation is valid at low energy. Let us justify this assumption in the next 
section. 



IV. CLASSICAL STABILITY I - RECOVERY OF FRIEDMANN EQUATION 

We can also show the essential, and somewhat surprising result, that the standard Fricdmann equation is recovered 
at low energy, starting from the action (5). There are higher derivative corrections to the Einstein equation due to 
the i?-dependence of and aF^. 

Since there arc higher-derivative terms, the stability of the system is a non-trivial question. In order to see the 
non-triviality, let us consider the Klein-Gordon equation (□ — M'^)ip = as a standard equation and add e{(fi)0^(fi/M^ 
to the right hand side. One might expect that the standard equation should be recovered whenever e — » 0. Actually, 
this is not true. The standard equation is certainly recovered in the e limit if e > 0. On the other hand, if e < 
then the system has a tachyonic degree and is unstable. Moreover, if ip crosses a root of e then the system experiences 
a singularity (□^<p diverges) and the low energy effective theory cannot be trusted unless a miracle cancellation occurs. 

In our system, from the estimate of Lkin in (17) with (7), Lkin is much smaller than aR^ (^ H^) at low energy 
if a 7^ 0. Hence, aR^ should dominate the higher derivative corrections and control the stability. For the theory 
roughly speaking, the parameter a plays the role of e above (including the sign) and it is known that the 
low energy dynamics is stable if and only if q > [34] . Here, stability means that as the universe expands, the system 
keeps away from unphysical spurious solutions and approaches the standard low energy evolution asymptotically. If 
we did not include the term aR^ then, as we shall see below, Lkin would make the quantity corresponding to e above 
to be negative essentially because f{R) is in the denominator. Therefore our system should be stable and the standard 
Priedmann equation should be recovered at low energy if and only if a > and (7) are satisfied. 

The equation of motion derived from the action (5) includes up to fourth order derivatives of the metric. However, 
since the Friedmann equation is a constraint equation, it does not include the highest order derivatives of the metric. 
This means that the generalized Priedmann equation describing a homogeneous, isotropic universe in our model 
includes derivatives of the metric only up to third order. Hence, it should be of the following form. 

e{t) ( H\ K?{pother + P4>fi + ^P4>)-^H^ .ION 

-Jfdt \H2j= ' (18) 

where the dimensionless coefficient e{t) is written in terms of {H/H"^, kH, k^tt), Pother is the energy density of fields 

other than (J), p^ Q = 2KLkin.K ~ Lkin -I- V is a part of the energy density of (j) which would be obtained by neglecting 
the i?-dependence of Lkin, ^Ptfi is corrections to p^^ which depends on H and H but not on H. We shall see the 
detailed form of the generalized Priedmann equation soon. 

In this form of the generalized Friedmann equation, e and Ap^ characterize corrections to the standard Friedmann 
equation. Indeed, if e = and if Ap^ = then this equation reduces to the standard Friedmann equation = 
ti^iPother + Pc/),o)- Howcvcr, this does NOT necessarily guarantees the recovery of the standard Friedmann equation in 
the limit e — > 0, Ap^ — > 0. As in the above example of a Klein-Gordon equation with a higher derivative correction, 
what makes the limit rather subtle is the sign of the coefficient e of the highest order derivative term. As we shall 
see soon, if e approaches zero from the positive side then the standard Friedmann equation is indeed recovered in the 
limit. This situation is somehow similar to having an infinitely massive extra degree of freedom. On the other hand, if 
e approaches zero from the negative side then the situation is similar to having an infinitely unstable tachyonic degree 
and, thus, the system is completely unstable. Moreover, if e reaches zero in finite time then the system encounters a 
singularity where dt{H/H^) diverges. Hence, both the stability and the recovery of the standard Friedmann equation 
require e to be positive. 

In our model there are two contributions to e: one from Lkin and the other from aR'^. We shall sec below that the 
first contribution is negative. This implies that the system would be unstable in the absence of the term aR^ and 
that we really need this term (or other ordinary higher derivative terms). On the other hand, we shall see below that 
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the sign of the second contribution is the same as the sign of a. Moreover, as we have already seen in the previous 
section, Lfei„ is much smaller than iJ^ and, thus, than aB^ at low energy. This means that the contribution from 
the olB? term determines the sign of e at low energy and that e is positive if and only if a is positive. Therefore, if 
and only if a is positive, the system is stable and the standard Friedmann equation is recovered at low energy. In the 
following, we shall show this statement more explicitly. 

For the stability analysis and numerical works, it is sometimes more convenient to work with a set of first order 
differential equations than a set of higher order differential equations. In the end of this section, we shall use the set 
of first-order equations to show the recovery of the standard Friedmann equation analytically. Some numerical results 
are shown in Sec. V. 

A set of first order differential equations is obtained in Appendix A for a general action of the form (Al) with 
(A24). Restricting it to our model action (5), we obtain 



7i- = -3il7r-y' ((/)), 

H = H^n, 

pother ~ — 3H{pother + Pother), 



(19) 



where 



H 

^ = 2{LRy - 2HLr, 
Lr = 2aR + Lkin,R, 



(20) 



and 



4a • 



(21) 



3kV 3(2g - 1) 



H 



+ 



[i2n{n + 2)H^f" - /'] {k\^)^ 



6{3H7T + V')f^,^ _ («4^2)-^ ^ 72aK'H\l - 2fl){n + 2), 



2n 



H 

1 



(Pother + Pother) + f ' (k TT ) ^"-^ 



3{2q l)nf^^,^,^_j^ ^ 72aK^H^n{n + 2), 



qn^ 



R=6H^{n + 2). 

Here, Pother and Pother are energy density and pressure of fields and/or matter other than (p. Besides the above 
dynamical equations, there is a constraint equation of the form 



H 



= :Fo + - maK^H'^{Q. + 2f, 



(22) 



where 



_ K'^iPother + P<P,o) " 3if^ 
-^0 = Jj2 ' 

0^,0 = 2KLkir,,K - Lkir, + V= " (kV^)^?^ + V. 



(23) 



Note that the standard Friedmann equation and the standard dynamical equation are = and Qq = Q, respectively. 
When the system is analyzed numerically, we can use the constraint equation to set the initial condition for y> and 
also to check numerical accuracy. 
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It is also possible to use the constraint equation to eliminate from the set of first-order equations: in the forth 
equation in (19) is rewritten as 



= jFo + AJT, (24) 



where 



Ajr = M?9_i) ^^20(0 + 2)H^f" - (Q + 1)/'] {k\^)^- 



H 



K^TT ■ (kV^)"^ - l08aK^H^n{n + 2). (25) 



Hence, the fourth equation in (19) is reduced to (18) with Ap^ = H'^ AT . 

Now let us show the recovery of the standard Friedmann equation at low energy by using the above set of first-order 
equations. 

From the inequalities in (17), both Lkin and oB^ (~ H^) arc small compared to the Einstein- Hilbert term R/2k^ 
(~ H^/k"^) at low energy (kH <C 1). Thus, even without any calculations we can conclude that 

e(i)-0, -^(p^,o-^)-0, ^L^=A^-0, (26) 

in the low energy limit kH — > 0. Of course, it is straightforward to confirm this by using the explicit expressions 
presented here. More specifically, the second term in the expression (21) of e{t) is small compared to the first term at 
low energy {kH <C 1): 

K-^f"{R) ■ (kV^)^ ~ (c2)55^(K2//2)2(m-l)-^ ^ (27) 

where we have used the condition (7). Note that the second term in (21) is negative if f"{R) is positive. Therefore, 
if and only if a is positive and the condition (7) is satisfied, e{t) is positive at low energy. In the next paragraph we 
shall see that the positivity of e{t) and, thus, the positivity of a and the condition (7) are essential for the recovery 

of the standard Friedmann equation. 

Having the behavior (26) and the positivity of e(i), it is easy to see that 

[K\pother + V)-iH^\/H'' (i^oo), (28) 

provided that H > and H/H^ < 0. This is equivalent to JFq ^ (t ^ oo) since the kinetic part in p^^ in (23) goes 
to zero much faster than the potential V. Actually, if the right hand side of (18) is positive (or negative) then H/H^ 
will increase (or decrease, respectively). Hence, H will decrease less (or more, respectively) rapidly. Moreover, pother 
decreases more (or less, respectively) rapidly if it satisfies the weak energy condition. Thus, as a result, the right hand 
side of (18) will decrease (or increase, respectively). The decrease (or increase, respectively) should continue until 
the right hand side of (18) becomes almost zero, and we expect the behavior (28), which is nothing but the recovery 
of the standard Friedmann equation at low energy. This behavior, namely the recovery oi the standard Friedmann 
equation, can be confirmed numerically. See Figs. 10-13 in the next section for numerical confirmation of the recovery 
of the standard Friedmann equation {Tq 0) and the standard dynamical equation {Qo ^0). In the low energy 
Friedmann equation, k^V plays the role of the cosmological constant Ag//. 



V. NUMERICAL RESULTS 



It is straightforward to integrate the set of first order equations (19) numerically. In this section we show some 
results. In the following we set 

q=l, m = 2, C = 1, a = 1; Pother = 0, Pother = 0, (29) 

for simplicity. In order to set initial values we need four independent conditions at i = since we have five variables 
and one constraint equation. In the following we show numerical results for four different sets of initial values. We 
shall see that the system exhibits an attractor behavior. Because of the attractor behavior, in figures showing results 
of long numerical integration, plots for those four different initial values are degenerate. 

First, let us confirm the behavior B = k^Htt/c —(3 — H / H'^)~^ stated just after (11). Note that the feedback 
mechanism and the stability analysis are based on this important behavior. Fig. 1 shows that —(3 — H/H^)B indeed 
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converges to 1 rather quickly. Thus, the evolution of tt is determined by this attractor behavior indecently of the 
initial condition. Figs 2 and 3 show the evolution of tt. 
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FIG. 1. The behavior of B = k^Htt/c for four different sets of initial values. In the figure, a is the scale factor and 
oo = a{t = 0). The initial values of <f> and H / arc ck(0 — 4>o)\t=Q = 0.01 and H /H'^\t={) — —0.02. Two more independent 
conditions at t = are given by setting J^o and Qo in four different ways as {J^o,Qo)t=o = (0.1,0.1), (—0.1,0.1), (0.1,-0.01), 
(—0.1, —0.01). This figure indeed shows that —(3 — H/H'^)B quickly converges to 1. 
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FIG. 2. The evolution of tt for four different sets of initial values. In the figure, a is the scale factor and oo = a{t = 0). The 
initial values of and H/H^ are ck((/> — (/>o)|t=o = 0.01 and H/H^\t=o = -0.02. Two more independent conditions at t = are 
given by setting J^o and Go in four different ways as {J^o, Qo)t=o = (0.1, 0.1), (-0.1, 0.1), (0.1, -0.01), (-0.1, -0.01). 
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FIG. 3. The asymptotic evolution of tt for four different sets of initial values. In the figure, a is the scale factor and 

ao = a{t = 0). The initial values of 4> and H/H^ are ck{4> — (f>o)\t=o = 0.01 and H/H'^\t=o ~ —0.02. Two more independent 
conditions at t = are given by setting J^o and Qo in four different ways as {J^o,Go)t=o = (0.1,0.1), (—0.1,0.1), (0.1, —0.01), 
(—0.1, —0.01). Four lines are degenerate because of the attractor behavior. 



Second, let us sec that the feedback mechanism indeed works and the effective cosmological constant stalls near 
zero. Figs. 4 and 5 show that. It is worth while stressing again that V{(j)) includes the ground state energies of all 
fields and that adding finite terms to the potential does never spoil the feedback mechanism as far as the minimum 
of V{(j)) remains negative. Although the additional finite terms change the form of V^((/>) but the behavior of T^((/>) 
near its root is always characterized by two parameters c and 4>o (i.e. the first order Taylor expansion). The feedback 
mechanism works irrespective of the values of c and ^o- 



0.01 



< 0.005 



(Fq, Go)t=o- 



(0.1, 0.1 ) 
(-0.1, 0.1 ) 
(0.1,-0.01) 
(-0.1,-0.01) 



I ' 1 

20 40 

ln(a/ao) 

FIG. 4. The evolution of K^Aeff = k^V{4>) — CK{(f> — (po) for four different sets of initial values. In the figure, a is the scale 
factor and ao = a{t = 0). The initial values of (f) and H/H^ arc CK{(f> — (f>o)\t=o = 0.01 and H /H^\t=o = —0.02. Two more 
independent conditions at t = are given by setting J^q and Qq in four different ways as (.7^o,So)t=o = (0.1,0.1), (—0.1,0.1), 
(0.1,-0.01), (-0.1,-0.01). 
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FIG. 5. The asymptotic evolution of K^Aeff = K^V{(j)) = ck{4> — (f>o) for four different sets of initial values. In the figure, 
a is the scale factor and oo = a{t = 0). The initial values of (f> and H/H^ are ck{4> — 4>o)\t=o = 0.01 and H /H^\t=o — —0.02. 
Two more independent conditions at t = are given by setting J^o and Qo in four different ways as {!FoiQo)t=o = (0.1,0.1), 
(—0.1, 0.1), (0.1, —0.01), (—0.1, —0.01). Four lines are degenerate because of the attractor behavior. 

Third, because of the behavior of Agyj governed by the feedback mechanism, the Hubble parameter H also ap- 
proaches zero very slovi^ly. Figs 6 and 7 show this behavior of H . Figs 8 and 9 show that the dimensionless evolution 
rate H /H^ of H indeed approaches zero. Thus, asymptotically H does not change in cosmological time scale. These 
again confirm the statement that Ae/ / stalls near zero and that it does more slowly than matter or radiation. 
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FIG. 6. The evolution of the Hubble parameter H for four different sets of initial values. In the figure, a is the scale factor 

and ao = a{t = 0). The initial values of </> and H /H^ are c«:(0 — (^o)|t=o = 0.01 and H /H'^\t={) = —0.02. Two more independent 
conditions at t = are given by setting To and Qo in four different ways as {ToiQo)t=o = (0.1,0.1), (—0.1,0.1), (0.1,-0.01), 
(-0.1,-0.01). 
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FIG. 7. The asymptotic evolution of the Hubble parameter H for four different sets of initial values. 

scale factor and 0,0 — a{t — 0). The initial values of (p and H/H^ are ck{(I) — 0())|t=() = 0.01 and H /H^\u 
independent conditions at t = are given by setting To and C/o in four different ways as {To,Qo)t=o = 
(0.1, —0.01), (—0.1, —0.01). Four lines are degenerate because of the attractor behavior. 
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FIG. 8. The evolution of the dimensionless evolution rate H /H^ of H for four different sets of initial values. In the figure, 
a is the scale factor and ao — a{t — 0). The initial values of and H/H'^ are ck{4> — <^o)li=o = 0.01 and H/H'^\t=o = —0.02. 
Two more independent conditions at t = are given by setting To and Go in four different ways as {T'o,Qo)t=o = (0.1,0.1), 
(-0.1,0.1), (0.1,-0.01), (-0.1,-0.01). 
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FIG. 9. The asymptotic evolution of the dimensionless evolution rate H / of H for four different sets of initial values. 
In the figure, a is the scale factor and ao = a{t = 0). The initial values of <f) and H/H^ are cK{cj} — 0o)|t=o ~ 0.01 and 
HIH\=q = -0.02. Two more independent conditions at t = are given by setting J^o and Qo in four different ways as 
i^o,Go)t=o = (0.1,0.1), (—0.1,0.1), (0.1,-0.01), (-0.1,-0.01). Four lines are degenerate because of the attractor behavior. 



Fourth, we can also see that the standard Friedmann equation is recovered (J-q 0) and that the standard 
dynamical equation is also recovered (Qq — > 0). What actually happens is that the standard dynamical equation is 
recovered earlier than the standard Friedmann equation and that the latter is gradually recovered after that. Figs. 10, 
11 and 12 show the early stage of the recovery process. Fig. 13 shows the late stage the recovery of the standard 
Friedmann equation. 
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FIG. 10. The early stage of the recovery of the standard Friedmann equation {J^o 0) for four different sets of initial 
values. In the figure, a is the scale factor and ao = a{t = 0). The initial values of (j) and H/H^ are ck{4> — <l>o)\t=o = 0.01 
and H/H^\t=o = —0.02. Two more independent conditions at t = are given by setting J'o and Go in four different ways as 
(:Fo,eo)t=o = (0.1,0.1), (-0.1,0.1), (0.1,-0.01), (-0.1,-0.01). 
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0) for four different sets of initial values. In the figure, a is the 



scale factor and ao = a{t = 0). The initial values of and H/H arc CK{(j) — = 0.01 and H /H |t=o = —0.02. Two more 

independent conditions at i = are given by setting J^o and Qo in four different ways as {J^o,Go)t=o = (0.1,0.1), (—0.1,0.1), 
(0.1,-0.01), (-0.1,-0.01). 
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FIG. 12. The early stage of the recovery process of the standard Friedmann and dynamical equations {J^o 0, Go 0) for 
four different sets of initial values. In the figure, a is the scale factor and ao = a{t = 0). The initial values of 4> and H/H^ are 



CK{(f> — = 0.01 and H /H^\t=o ~ —0.02. Two more independent conditions at t 

four different ways as {To,go)t=o = (0.1,0.1), (-0.1,0.1), (0.1, -0.01), (-0.1, -0.01). 
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FIG. 13. The late stage of the recovery of the standard FYiedmann equation (J^o — > 0) for four different sets of initial values. 
In the figure, a is the scale factor and oo — a{t = 0). The initial values of (f> and H/H^ arc ck{4> — (^o)|t=() ~ 0.01 and 
H/H^\t=o = -0.02. Two more independent conditions at t = are given by setting J-'o and Qo in four different ways as 
i^o,Go)t=o = (0.1,0.1), (—0.1,0.1), (0.1,-0.01), (-0.1,-0.01). Four lines are degenerate because of the attractor behavior. 

Finally, we can check that the extra degree of freedom if goes to zero. Figs 14 and 15 confirm that behavior. 
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FIG. 14. The evolution of the extra degree ip for four different sets of initial values. In the figure, a is the scale factor and 

ao = a{t — 0). The initial values of (/> and H/H^ arc ck{(I> — 0o)|t=o ~ 0.01 and H/H'^\t=o — —0.02. Two more independent 
conditions at t = are given by setting J^o and Qo in four different ways as {J^o,Go)t=o = (0.1,0.1), (—0.1,0.1), (0.1, —0.01), 
(-0.1,-0.01). 
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FIG. 15. The asymptotic evolution of the extra degree (p for four different sets of initial values. In the figure, a is the scale 
factor and ao = a{t = 0). The initial values of 4> and H/H^ arc CK{(f> — 4>o)\t=o = 0.01 and H /H^\t=a — —0.02. Two more 
independent conditions at t = are given by setting To and Qo in four different ways as (To,Qo)t=o = (0.1,0.1), (—0.1,0.1), 
(0.1, —0.01), (—0.1, —0.01). Four lines are degenerate because of the attractor behavior. 

VI. REHEATING 

We have achieved the vanishing cosmological constant in a way that is stable under radiative corrections and that 
has self-consistent, stable dynamics. However, although the cosmological constant approaches zero, it docs so more 
slowly than matter or radiation so that without additional structure, the universe would be empty. It is not entirely 
clear that a dynamical model where this is not the case could be successful since it is this property that makes it 
possible for all fields other than (p to settle into their ground state before (f> stalls at zero curvature so that the zero 
curvature really corresponds to the vanishing cosmological constant. Moreover, the singular behavior of a kinetic 
term coefficient and, thus, the slow evolution of (f) are required for stability. This does imply, however, that should 
this mechanism be responsible for a low cosmological constant, reheating would be required to thermally populate 
the universe after the cosmological constant has decreased to a small value. Nonetheless, as we shall see below, this 
model can at the very least reduce fine-tuning by 60 orders of magnitude or provide a new mechanism for sampling 
possible cosmological constants and implementing the anthropic principle. 

The reheat process requires further speculation, and is a subject for future research. A couple of possibilities are: 

(I) Low-energy inflation. One can consider an extra scalar field x with mass ~ lO^^eV and a term like —Rx^. 
When R ~ m^, a phase transition would occur (as in hybrid inflation [35]) and the universe would be reheated 
up to temperature ~ TeV. This phase transition happens when the energy stored in % plus the energy stored in 
<i> yields a Hubble constant of approximately m-^. The energy in x '^'iH decrease during the phase transition; the 
energy after the phase transition must be very small. In this case the cosmological constant problem is reduced from 
{Mpi/lQ-^eVy ~ 10^2° to {TeV/lQ-^eVy ~ 10^°. For smaller m^, the reheat temperature would be lower and the 
tuning of the cosmological constant would presumably be smaller. 

(II) Energy inflow from extra dimensions. For example, in a non-elastic scattering of branes, a part of the kinetic 
energy due to the relative motion can be converted to radiation on our brane without changing the brane tension and 
the cosmological constant. For this to work, branes should be sufficiently flat and parallel. 

After reheating, the conventional standard cosmology with a vanishingly small cosmological constant can be restored. 
This is because the scalar field ^ is just frozen at low energy. Namely, 

\k'^4>\ ~ (kF)^"^ X \k^tt\^ (30) 

can be made arbitrarily small at low energy {kH <C 1) by considering a sufficiently large m, where k^tt is estimated 
by (12). In the scenario (I), both n and H evolve continuously through the reheating epoch. Hence we obtain 

\k^V\ ^ (kH)^""-^ , (31) 
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and the 4d cosmological constant does not overshoot zero within the present age of the universe if m is large enough 

to ensure that (K^i?^^^)^™- 2(2,-1) < KHtoday[K'^Rreh/4: - K'^{U{x+)\R=R^^h - f^(X-))]) where Rreh is the value of the 
Ricci scalar at reheating, x+ and X- are the expectation value of x before and after the phase transition, respectively. 

Hereafter, we set c = 0(1). 

On the other hand, in the scenario (II), while tt evolves continuously, H increases suddenly by energy flow from the 
extra dimension. Hence, H in (12) and (30) should be estimated before and after reheating, respectively. Therefore, 
we obtain from (30) the following estimate for the period soon after reheating. 

k\V\/V ~ {KHafter)^"'/{KHt„fore)^^\ (32) 

where we have used the Friedmann equation •iH^f.fore k'^V. Here, Hbefore and Hafter are the Hubble parameter 
before and after the energy inflow from extra dimensions, respectively. By using (11) it is easily shown that k^ttH + 
c/5 oc during the radiation-dominated epoch {H / H"^ = —2). Hence, k^ttH decays from ^ Hafter /Htef ore to 0(1) 
in the time scale Atreiax ~ {Hafter /Hbefore)^^^H~^^^^. Thus, if m is large enough to ensure that {KHafterY™'~^ < 

{KHbefore)'^^~^~ then V does not overshoot zero in the time scale Atreiax, provided that KHafter *C 1. After that, 
we can use the estimate (31) with H being the realtime value, and V does not overshoot zero if m is sufficiently large. 
Let us recall that a choice like f{R) ~ ex.p[— k~'^ R~^] corresponds to m — » oo. 

One might worry about symmetry restoration and phase transitions that occur after reheating. This is not a 
problem in both examples if m is sufficiently large. The reason is that Rreh/4: — '«^(C^(x+) h ^ U{x-)) (or k^V, 
respectively) is the cosmological constant at zero temperature since the temperature before the reheating is zero. The 
large m ensures that Atoday is still positive and small since (p continues to be almost frozen all the way down to the 
present epoch including the time when the symmetry is restored and during the time the phase transition takes place. 

Even if all else fails, although not our initial subjective, the existence of (j) can at the very least provide a natural 
framework in which to implement the anthropic principle. If wc assume eternal inflation, there would be many 
inflationary universes. In each universe, the cosmological constant is determined by how much has rolled when 
inflation ends. That in turn depends on the number of e-foldings that occurred before inflation stopped. In an eternal 
inflation scenario, different numbers of e-foldings would occur in different domains and therefore different (j) values, 
and hence different values of the cosmological constant would occur in different regions. 

Our model in general predicts = p<f,/ p<i) — —1 today because of (31) since pre- and post-reheating behavior re- 
quires a large m. In the scenario (I), \ V\ continuously decreases. In the scenario (II), \ V\ increases suddenly at reheating 
but is adjusted to the behavior (31) in the time scale Atreiax- Note that HtodayAtreiax < {Hufore/ HafterY^^ 1 
since Hl^^^^ > H^^fore — ^before/^ > ^today/^ — Hf^^ay This means that the behavior (31) and ~ -1 are 
realized before the present stage of the universe. 

VII. CLASSICAL STABILITY II - LINEARIZED GRAVITY 

In this section, restricting to the q = 1 case for simplicity, we show that linearized Einstein gravity in a Minkowski 
background is recovered at distances longer than the length scale /* = \/aK and at energies lower than The 
extension to general q (> 1/2) should be straightforward. 

A possible source of instability which would disturb weak gravity in a Minkowski backgromid could be an inho- 
mogeneous fluctuation of ^, which could possibly generate violent breakdown of linearized Einstein gravity. Quite 
surprisingly, this is not the case and the seemingly most dangerous part, the kinetic term, is not as dangerous as 
it looks. Essential to this conclusion is the constraint equation, which prevents (j) from fluctuating freely and forces 
the denominator and the numerator to fluctuate in a strongly correlated way so that the contributions of the kinetic 
term to the equation of motion are regular and much smaller than those of the aR^ term. Namely, the scalar field 
and metric cannot fluctuate independently. In fact, it is a well-known fact in the standard scalar fleld cosmology that 
there is only one physical degree of freedom out of scalar field and metric degrees for scalar- type perturbations. In 
the longitudinal gauge (see (40) or (Dll)), one of the constraint equations is the (Oi)-component of the perturbed 
Einstein equation, which is roughly of the form 

7r5,(<5</.)=.^-<5r(f'^^'-\ (33) 

where tt is the momentum conjugate to the homogeneous background (f) defined in (9), S4> is the inhomogeneous 
perturbation of 0, dT^"^'^'^^^ is the (Oi) -component of the stress-energy tensor of fields and matter other than </>. Note 
that homogeneous perturbations can be absorbed into the background and that, without loss of generality, we can 
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restrict our analysis to inhomogeneous perturbations, for which di{5cf)) 7^ 0. The right hand side of this equation is 
manifestly regular and we already know that tt behaves like (12). Thus, we conclude that 

6(j) ~ kH X (regular expression). (34) 

This means that is indeed suppressed at low energy (kH <^ 1). Similarly, the other constraint equation, namely 
the traceless part of the (ij)-components of the perturbed Einstein equation, leads to 

(50)' = {kH)^ X (regular expression). (35) 

This again means that inhomogeneous perturbation of (p is suppressed at low energy. Hence, from the constraint 
equations, we expect that the stress energy tensor of the scalar field should be small enough. More rigorous treatment 
of the constraint equations is included in the complete analysis shown in Appendix D. 

Since contributions of the singular- looking kinetic term are expected to be small enough, the only possibly important 
correction to Einstein gravity is again due to the aB? term. This tells us that the linearized gravity in our model 
should be similar to that in the theory R/2k^ + aR^ . Hence, we should be able to recover the linearized Einstein 
gravity in Minkowski background at distances longer than the length scale /* = \/aK and at energies lower than 
[36]. 

A complete analysis of the weak gravity in Minkowski background requires much more carefulness. One of the 
reasons is that the action for the scalar field is not manifestly well defined in Minkowski background since the 
denominator of the kinetic term vanishes at i? = 0. Moreover, it is not totally clear how to treat the perturbations of 
the denominator. This situation requires a sort of regularization for the field equation. In the end of the calculation 
we of course need to take the limit where the regularization is turned off. Unlike calculations in quantum field 
theories, we must not renormalize anything since we are dealing with classical dynamics. Namely, before turning off 
the regularization, we must not subtract anything from the regularized field equation and must treat it as it is. 

For the purpose of the rigorous treatment of weak gravity, we investigate in detail perturbations around the flat 
FRW background in the longitudinal gauge and take the kH +0 limit in the end of the calculation. This strategy 
makes it possible to analyze linear perturbations around Minkowski background, on which the action of the scalar 
field perturbation is not apparently well-defined. For simplicity, we assume that the FRW background is driven by the 
scalar field cf). Of course, we include a general matter stress-energy tensor into perturbations and carefully investigate 
how they couple to gravity. A critical point for this analysis is that constraint equations ((0«)-components and the 
traceless part of (ij )-components of Einstein equation) prevent </> from fluctuating freely, and it turns out that the (p 
perturbation is of sufficiently high order in the kH expansion that the perturbative analysis is under control. 

The result of the analysis is simple. The linearized gravity in our model is similar to that in the theory R/2k^ + aR'^ 
in the sense that the only difference is the existence of an extra scalar-type massless mode in our model. The extra 
massless mode is decoupled from the matter stress-energy tensor at the linearized level. Hence, as far as we are 
concerned with classical, linear perturbations generated by matter sources, these two theories give exactly the same 
prediction. Therefore;, in our model the linearized Einstein gravity in Minkowski background is recovered at distances 
longer than the length scale U = -Jolk and at energies lower than . 

Let us begin by deflning the metric perturbation bg^j.^ by 

9ixi.=9f}^^9ixv, (36) 

where 

g^^^dx^dx" = -dt^ + a{tfSijdx'dx^ (37) 

{i = 1,2,3). We shall investigate the system of coupled equations for the scalar field perturbation S(f), the metric 
perturbation Sg^^i, and the stress energy tensor contributions from other fields and matter, taking into account all 
higher derivative corrections. In order to investigate linearized gravity in Minkowski background, in the end of 
calculations, we shall take the limit 

kH +0, (38) 

keeping H / H"^^ {n — 1, 2, 3) finite, where H — dta/a. In the following, we shall denote dt by an overdot. 

In order to take the full advantage of the symmetry of the background spacetime (with iJ ^ 0), we expand the 
scalar field perturbation, the metric perturbation and the stress energy tensor contributions from other matter and 
fields by harmonics on the 3-plane with each expansion coefficient being a function of the time t only. Now it is 
well-known and easily shown that there are three distinct types of perturbations and that each type is decoupled 
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from the others in the Hncarizcd level. They are called scalar-, vector- and tensor-type perturbations and each of 
them is an irreducible representation of the symmetry of the background FRW spacetime: spin-0, spin-1 and spin-2 
representations, respectively. Since scalar-, vector-, and tensor- type perturbations are decoupled from each other in 
the linearized level, we can analyze perturbations of each type separately. 

The vector- and tensor-type perturbations are decoupled from scalar-type perturbations and, thus, from the pertur- 
bation of the scalar field (j). Moreover, the linear perturbation of the Ricci scalar vanishes for vector- and tensor-type 
perturbations. Thus, the kH +Q limit of perturbation equations is manifestly well-defined. Actually, in the 
kH +0 limit, the linear perturbation of the stress energy tensor T^^^ of (p vanishes. Therefore, we obtain 

6G^. = «'t;*''^^ (39) 

This is nothing but the linearized Einstein equation. In other words, possible corrections to the linearized Einstein 
equation appear only for the scalar-type perturbations. 

So, let us concentrate on the scalar-type perturbations. As is well-known, in linear perturbations we have not 
only physical degrees of freedom but also gauge degrees of freedom so we need to fix gauge. (See subsection D 1 of 
Appendix D for the formula of infinitesimal gauge transformation and a gauge choice.) In the longitudinal gauge, the 
perturbed metric is written as 

g^^dxi'dx" = -(1 + 2^Y)dt^ + 2^Y)a^5ijdx'dx\ (40) 

where $ and are functions of t, and Y is the scalar harmonics on the 3-plane with 3-momentum fcj. Here, we omit 

to write ki and the integration over all possible values of ki. 

As already stated, we shall analyze the system of coupled equations for the scalar field perturbation 5(j), the metric 
perturbation 5gij,v and the stress energy tensor contributions from other fields and matter in the FRW background, 
taking into account all higher derivative corrections. In the end of calculations, we shall take the limit kH +0, 
keeping i9"i//7J"+^ (n = 1, 2, 3) finite. For this purpose, we introduce a small dimensionless parameter e so that 

kH = 0(e), 



0(e°). (41) 

and expand all quantities and equations. We shall keep the expanded equations up to the order O(e^) to obtain 
equations governing the quantities of order 0(e°). In the end of the calculation, we shall take the limit e ^ 0. 

Complete analysis of the scalar- type perturbations is given in subsection D2 of Appendix D. The equations 
governing the 0(e°) order metric perturbations are obtained after going through the procedure described above and 
eliminating the scalar field perturbation. The result is summarized as the following equations for linearized gravity 
in Minkowski background. 



«;2fc2*+ + (too + 2fcV(ii)) = 0, 



□ 



{l-12aK^a)^_+4:aK^ (3<i/+ + k^'i/+^ -2k^^ll) =0, (42) 



where t's are scalar-type perturbations of the stress energy tensor contributions from other fields and matter defined 
by (D20). Here, we have redefined the spatial coordinates as a^x^ x^ and 

□ = - A:^. (43) 

Before taking the kH limit, oq is defined as the value of the scale factor at the time when a weak gravity 
experiment is performed, assuming that the duration of the experiment is much shorter than the cosmological time 
scale H^^. After taking the kH limit, the scale factor is a constant and, thus, the above rescaling of the 
spatial coordinates makes perfect sense at any time. The perturbed metric, after the kH +0 limit and the above 
redefinition of the spatial coordinates, is 

ds^ = -(1 + 2<i>Y)df + {1- 2^Y)d^jdx'dx^, (44) 

and \l/± = * ± The conservation equation is 

Too + fc^T(i)0 = 0, 

T(L)0 + ^I^^T{LL) - T(y) = 0. (45) 
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Now it is time to show that the hnearized Einstein gravity in Minkowski background is recovered at distances longer 
than the length scale U = \/aK and at energies lower than l^^. (Note that we have assumed that a is positive.) For 
this purpose and to make the arguments qualitative, let us compare the result (42) for scalar-type perturbations in 
our model with the corresponding equations in the higher-curvature theory whose gravitational action is 

Ihd = J d^x^g + olR^ ■ (46) 

For scalar perturbations given by (44) and (D20), linearized gravity equations in the theory (46) are 

+ (too + 2k'^T(LL)) = 0, 

(1 - 12aK^n) + Aan^ + A;^*+) - 2k^t(ll) = 0. (47) 

Therefore, for linearized gravity the only difference between our model and the higher-curvature theory (46) is that 
there is an extra massless mode for in our model. (Notice the extra □ in the second equation in (42).) However, 
this zero mode does not couple to the matter stress energy T°*^'^^ directly since the extra □ in the second equation 
in (42) is applied to the whole expression. Hence, as far as we are concerned with linear perturbations generated 
by matter sources, these two theories give exactly the same prediction. Thus, it is concluded that the scalar-type 
linearized gravity in Minkowski background in our model is effectively described by the theory (46). 

In the theory (46), the linearized Einstein gravity is recovered at distances longer than the length scale = ^/aK 
and at energies lower than To see this, let us quote the equations governing scalar perturbations in Einstein 
gravity: 

K''kH+ + (too + 2fc2r(ii)) = 0, 

- 2k\ll^ = 0. (48) 

It is easy to see that these equations in Einstein theory are recovered from the corresponding equations (47) in 
the theory (46) if q:k^4'±/4'± and an^k'^ are sufficiently small. For the recovery, it is fairy important that we had 
assumed that a is positive. Actually, the higher derivative terms in the second equation in (47) do not introduce 
extra instabilities if and only if a is non-negative. Note that the stability of homogeneous, isotropic evolution of the 
universe has also required a > 0. 



VIII. QUANTUM MECHANICAL STABILITY 

Now let us think about quantum mechanical stabilities. 

The stability under radiative corrections is an important feature of our model. Radiative corrections will produce 
additional regular terms in the action, but the field will stall whether or not these are present. Adding finite terms to 
the potential part does not change anything since the curvature R feels the corrected potential via the gravity equation. 
As for the kinetic part, we would like to stress again that the singular behavior of the kinetic term coefficient is imposed 
only on the most singular- looking term among many possible terms in the kinetic part and, thus, adding any kinetic 
terms which are less singular-looking at i? = does not change anything. Of course, adding a more singular-looking 
kinetic term just makes the condition more robust. The more singular a kinetic term looks, the more stable it is under 
radiative corrections. 

Since /(i?) is in the denominator of a kinetic term and vanishes at i? = 0, one might also worry about additional 
singular potential terms like 1/ f{R) being generated through radiative corrections. However, this does not happen. 
In order to see this, it is convenient to normalize quantum fluctuation of the scalar field </> around the homogeneous 
classical background. For simplicity, we restrict our discussion to the q = 1 case and work in the unit where k = 1. 
The kinetic term is expanded as 



2/ 



1 



{SR + 62R) + 



{SRy 



(^r - h"" + h'^^K) {<Po + V^<5</>c) (00 + ^/ToHc) 



(49) 



up to the quadratic order in perturbations, where a quantity with the subscript represents the value on the FRW 
background considered in Sec. IV, /i^i^ is metric perturbation, 6R and ^2-^ are the 0{h) and 0{h?) parts of the Ricci 
scalar perturbation, respectively, and 
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is the normalized quantum fluctuation of the scalar field. Hence, the contributions of the kinetic term to the full 
quadratic action multiplied by —2 are 



JO -tU) 4 \ /o -n-O / -"'0 \J / ^ 



4 I 


fo 


Rofo 


Rq 


fo 


Rq 



-"O \//o /O -n-O -n-o Jo -Ko V -"-0 / 



Actually, all except for the first term vanish in the low energy limit since 

RoocH^ ^ 0, 

^0 rr H 



i?0 



oc H^""-^ (52) 



in the low energy limit if — > and 



- Roh(^)' -Rlfoi-,] (53) 



u,/vji pi U./'-'l P 

\J/o VJ 



This means that 6(j)c asymptotically approaches a canonically normalized fluctuation in the low energy limit. In terms 
of 5(j)c, each term in the potential part can include a positive power of \/Jo but not negative power. Therefore, loop 
contributions of 5<pc to the radiatively corrected potential part can include positive powers of vTo but not negative 
power. 

More rigorous treatment requires careful consideration of constraint equations among the scalar field perturbation 
and the metric perturbation. Namely, we need to introduce a variable analogous to the Mukhanov variable [37] 
in the standard field theory cosmology. Hence, rigorous treatment seems miicli more complicated than the above. 
Nonetheless, the above argument is convincing enough and we do not expect terms with negative power of y/Jo to 
appear in the potential part of the proper perturbation variable. 

In the above, it has been shown that the feedback mechanism is stable under radiative corrections and thus a zero 
or small cosmological constant is protected against radiative corrections. In the following, we show that quantum 
fluctuation at low energy is so small that the effective cosmological constant at low energy does not overshoot zero 
even quantum mechanically. 

In the K^i? — > limit the scalar field becomes completely weekly coupled. Hence, one might think that the scalar field 
could overshoot zero curvature by quantum fluctuation. This does not happen, as far as the (effective) cosmological 
constant is a substantial component of the energy density of the universe. The reason is as follows. (For simplicity 
we consider the q = 1 case only, but it is easy to generalize it to a general q > 1/2. We shall again work in the unit 
where k = 1, and assume that the stability condition m > 3/2 is satisfied.) 

The fluctuation 6(j) of the scalar field (p around a homogeneous background does not have a canonically normalized 
kinetic term. As we showed in the above, the normalized fiuctuation 6(pc = S(f)/\/Jo has a canonically normalized 
kinetic term plus additional terms which vanish in the low energy limit H ^ 0. Hence, by dimensionality, amplitude 
of the quantum fluctuation 6(j)c should be \6(j)c\ ~ H. Hence, the quantum fluctuation of the cosmological constant is 
estimated as 

\6Aeff\ ~ \cS4>\ oc ij2m+l^ (54) 

where we have used the behavior /o a H'^"^ near H = 0. On the other hand, if the (effective) cosmological constant 
is a substantial component of the cosmological energy density then 

Ae// ~ H\ (55) 
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Note that under the stabiHty condition m > 3/2 (see (7) with q = 1), the cosmological constant asymptoticaUy 
dominates the cosmological energy density. Since 2m + 1 > 2 from the stability condition, the estimates (54) and (55) 
imply that 

|<5Ae//|«Ae/J (56) 

at low energy {H <^ 1). This means that Ae// does not jump to a negative value by quantum fluctuation. 



IX. SUMMARY 



In the present paper, we have investigated gravity in the recently proposed dynamical approach to the cosmological 
constant. We have shown that (i) the effective cosmological constant decreases in time and asymptotically approaches 
zero from above; (ii) the evolution of a homogeneous, isotropic universe is described by the standard Friedmann 
equation at low energy; that (iii) classical, linearized gravity in Minkowski background is described by Einstein 
gravity at distances longer than = ^/an and at energies lower than where k is the Planck length and a is a 
dimensionless, positive parameter of the model; that (iv) the mechanism is stable under radiative corrections and thus 
a zero or small cosmological constant is protected against radiative corrections; and that (v) quantum fluctuation at 
low energy is so small that the effective cosmological constant at low energy does not overshoot zero even quantum 
mechanically. 

One of the most disturbing difficulties in thinking about the cosmological constant is that it is not protected against 
radiative corrections, which usually generate enormous vacuum energies compared to what wo observe. Because of 
the above properties (i)-(v), the feedback mechanism can be considered as a dynamical way to protect a zero or small 
cosmological constant against radiative corrections. Hence, although the feedback mechanism by itself does not solve 
the cosmological constant problem, it can help solving the problem. 
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APPENDIX A: DERIVATION OF EQUATIONS OF MOTION 



In this appendix we consider a more general action of the form 

R 



-9 



2k2 



L{R,X,K,cP) 



(Al) 



where X = W^Rij,^, K = -K^d^'ctydf^cj), and R and R are the Ricci scalar and the Ricci tensor of the metric Qj_tj/- 
(See Appendix B for an alternative description.) Our sign convention for the metric is ( — h ++). 



1. Variational formulas 



In this subsection we derive some general formulas for variations of geometrical objects in £)-dimension. 
Let us decompose the metric g^v into the background g^S and perturbation Sg/j,,^: 

9f.. = 9^^} + 5gi,u. (A2) 

The indices of any linear-order quantities are lowered by g^^} and raised by the inverse g^^^i^" of g^^}. For example, 

5gi^- = ^^"^''^^^"^'^"^a. (A3) 
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Hence, 

gi"" = g^^^"" - Sg"" (A4) 

up to the linear order. 

Next, we can easily expand \/— 5 &iid the Christoffel symbol T^,^ as follows, where g is the determinant of g^v 

- (1 + Isg^ , 

r^. = r(°);, + 5r^, (A5) 

up to the linear order, where Sg = g^'^^^'^Sg^^, r^°'^^ is the Christoffel symbol for the background metric gjS), and 

= Ig^^^^'iSga^;,^ + Sga.;^. - (A6) 

Here, a semicolon denotes the covariant derivative compatible with the background metric g^S} ■ 
Thirdly, the Ricci tensor is expanded as follows. 

up to the linear order, where is the Ricci tensor for the background metric gj^J , and 



Therefore, we obtain 



where 



= ^(^fl^;>.p + ^9^;f.p - ^9p^:fp - ^a-.v^u)- (A8) 



R = + 6R, 

X = + SX, (A9) 



6R = -R^^^'^'^Sg^. + {Sg^''., - 69'^;^, 

5X = 2i?(°)''''^ii^^ - 2i?(°)wij(0)^^^^^_ (Aio) 



2. Equations of motion 

By using the variational formulas presented in the previous subsection, the variation of the action is calculated 



as 



where 



(All) 



+2(L.xi?^''):^ + 2(L,xR''nfp - 2{^,xRnfp - 2{LxRn;p<rg'"' - ^L^xR^^'R;, 
E,i. = 2K\L,Kd^^y^' + L^. 

Hence we obtain the following set of equations of motion. 

where T'^/^er stress energy tensor of other fields whose action should be added to the action (Al). 



(A12) 



(A13) 
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Let us consider a homogeneous = (p{t) in the D = A flat FRW background (8). With this ansatz, the stress energy 
tensor is 



rplX 



( -P<t> Q 

P0 

P0 

V 



(A14) 



where 



P<t'+P<t' = 2K;*L,if(^^ + 2 
and 



L K + 2(3if2 + 2ij-)L X -2if(Lfl + 6fl'%x)- +4ij L « + 6(if2 + £r)L x , (A15) 



R = Q{H + 2H^), 

X = 12(ij2 + 3i/2^ + 3i/4^ 

Here, a dot represents the derivative with respect to the proper time t. The equations of motion are 

^2 

-ff = {P<j> + Pother) , 
-ff = - Y [('°'/' + P'^') + (Pother + Pother)] , 



= TT + SHn - L,, 
and the equation of motion for other fields, where 



TT = 2k*L 



and Pother{t) and Pother (t) are energy density and pressure of other fields: 

/ -Pother 

Pother 
Pother 

\ Pother 



T ^ = 

other V 



The conservation of T^^^^^ ^ is expressed as 

Pother + 'iH{pother + Pother) = 0. 

Note that (A17), (A18) and (A19) are not independent. Actually, there is an identity 

P0 + 3H{p4, + p^) + ^E^ = 0, 
which can be checked explicitly and corresponds to the conservation equation of T^^^. 



(A16) 



(A17) 

(A18) 
(A19) 



(A20) 



(A21) 



(A22) 



(A23) 



3. First order equations 

Hereafter in this appendix, we assume that L is of the form 

L = ^F{R, X)K^ - Vi^) + G{R, X), 
where F and G are functions of R and X, and g is a constant. In this case, p^ and are expressed as 



(A24) 
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q 



where 



+6K^/f2f2 + 2H'^{2Q. + 3)G,x] + kV, (A25) 

= 2[L,i?. + 2(3i?2 + 2ij)L,x]' - 2iJ(L,7j + 6ff^L,x). (A26) 
Hence, the equations of motion (A18), (A19) and (A22) are rewritten as the following set of first-order equations: 

^= fn ■ (K''7r2)2(2g-i)~5^ 
7T = -SHtt -V'{(l)), 

H = H^n, 

pother — -SH{pother+ Pother), (A27) 

where 

e = Uk^H^ [3G,flfl + 12H^{2Q + 3)G,rx + 12H^{2a + 3fG,xx + 2G x] 
6(2g-l)iJ2 



qK^ 



^3(2g_2) ^ 2(4fi2 + 6fi - 3)/fV,x 



^0 = 212 



+12Q [(O + 2)H^f^RR + 2(40^ + 130 + 12)H^f^Rx + 8(20 + 3)(02 + 3f2 + 3)H^f^xx] } (k^tt^)^ 
-6k2 {-G,ij + 2(4Q2 + 60 - 3)H^G,x 

+120 [(O + 2)H'^G,RR + 2(40^ + 130 + 12)if^G,ijx + 8(20 + 3)(02 + 30 + 3)iI^G,xx] } , 
1 



2^2 



K (Pother + Pother) + f ' [k TT )2<!-i 



AG EE _'^^:^!L^E + 2H^2n + 3)/,x] («V')^ + 6n'n [G,« + 2ir2(2J7 + 3)G,x] , 
i? = 6F2(0 + 2), 

X = 12H'^{n^ +3n + 3). (A28) 
The constraint equation (A17) is now written as 

-^=^o + Aa--^, (A29) 

where 

•^0 = Jp ' 

p^,o = 2KLkir.,K - Lkin + V= • (kV)^ + F. (A30) 



25 



When the system is analyzed numerically, we can use the constraint equation to set the initial condition for and 
also to check numerical accuracy. 

It is also possible to use the constraint equation to eliminate ^ from the set of first-order equations: T in the forth 
equation in (A27) is rewritten as 

= jTo + AJT, (A31) 

where 

AJ^ = _&Lp:l [j^^ + 2H\2n + 3)/,x] K^TT ■ (ac^tt^)-*^ 
^3(2Q_a) ^ ^^^^^ ^ 2(202 ^ _ ^)H^f ^ 

+12Q [(O + 2)H'^f^RR + 2(4f2^ + 13f2 + 12)H^Irx + 8(2f2 + 3){n'^ + 30 + 3)H^f,xx] } iK\'^)^ 
-6k^ {-(O + + 2(20^ + 30 - 3)iI^G,x 

+120 [(O + 2)H^G^RR + 2(40^ + 130 + 12)H'^G^rx + 8(20 + 3)(0^ + 30 + i)H^G^xx\ } - -tt^-- (A32) 



APPENDIX B: EQUIVALENT ACTION 

In this appendix let us consider an action of the form 



(Bl) 



where R is the Ricci scalar of the metric gf^^, (p is a, scalar field and K = —n'^g^'^df^cpd^cj). We shall show that this 
action is equivalent to the following action. 



\ ) 2k^ 

where the new field V' and the function -ftT, 0) are defined by 

eT«'^ = 2k^-¥r{R, K,(I))^R = R{ip, K, <j)), 

the metric is defined by 



(B2) 



(B3) 



(B4) 



and R is the Ricci scalar for the new metric 5^1/. Here, k is an arbitrary positive constant, = {dF/dR)K,(p, 
K = —K^g^^^d^4>d„(j), gi^" = and constants a, (3 and 7 are 



D 



a = 



7: 



^{D-l){D-2y 
2 

^{D-l){D-2y 



D-2 



D - 1 



(B5) 



In (B3) the function R{^, K, 0) is defined by solving the left equation with respect to R and, thus, we have implicitly 
assumed that d^F/dR? ^ 0. The action (B2) actually describes the Einstein gravity plus two scalar fields (j) and V'- 
Equations of motion derived from the two actions are the same and are 
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yZ^^^-e-i^oVp^a^^j + _e-"«'/'F^ = 0, (B6) 



where it is understood that R = e'^'^'I'k, </>) and K = e^^'^K are substituted into R, F, = {d¥/dK)R^^ and 
F^ = {dF/d(j))n.K- The equation of motion for can be derived from the second equation and the divergence of the 
first equation, and is 



gg^^d^tp] - a«;e-"'^'^F - 0K^e-^'^'^FKg^''di,(t>d^(l) + ^e'^'^'^R = 0. 



(B7) 



Again, it is understood that R = R{i}}, e^'^'^'i^, cj)) and K = e^'^'^ii' arc substituted into R, F and F^- 

Now let us show that the actions (Bl) and (B2) are classically equivalent. It is an extension of ref. [38], in which 

the scalar field has a usual canonical kinetic term. 
First, it is easy to show that 



-L(5 {^f^F{R, K, ct>)) = \Fg^'''5g^, + FrSR + FkSK + F 
/ g I 

= E^''Sg^^ + E5cP + X^.^, 



(B8) 



where 



E = 2K\FKd^,4>Y^ + F^, 

= FniSg'^-'., - Sg-',''') - FR.^Sg^'' + F^^Sg-", - 2k''F K<f>'^ 6 cj,, 



(B9) 



and a semicolon represents the covariant derivative compatible with the metric 9^,^. Equations of motion derived from 
61 ~ are E'^'^ — and E — 0. These equations, supplemented by suitable initial conditions, govern dynamics of 
the system. Evidently, unless F^ is independent of R, E^'^ includes up to the forth order derivatives of the metric 
gp^i,. In the following, we shall introduce an auxiliary field so that the resulting action includes only up to the second 
derivatives of fields. Moreover, it will be shown that after a conformal transformation, the gravitational part of the 
resulting action is of the form of the Einstein-Hilbert action. 
Second, let us perform a conformal transformation 



2a) 

g^v — 6 9tivi 



(BIO) 



where a; is a function to be determined below. The relation between the Ricci tensor i?^^ of the original metric g^i, and 
the Ricci tensor .R^y of the conformally transformed metric g^v can be found in textbooks in general relativity [39] . 



R^^i. = R-t,. -{D- 2)u),^uj,, + {D- 2)uj'Puj,pg^^ + {D - 2)a;.^, + uj-'P^g^^. 
Accordingly, the relation between Einstein tensors G^^ and G^^ for g^^, and g^j,, respectively, is 



Gij,v — Guv ~ {D — 2) 



+ {D-2){w.,^^-w'Ppg^^). 



(Bll) 



(B12) 



Thus, we obtain the following expression of Eij,,^. 



1 



-E^, = FflG^, - ~ Ffli?)5^, -{D- 2)Fr 



+{D - 2)Fr{uj,^, - uj'Ppg^^) - {Fr.^^ - F/.^g^,) 



1 



K^F 



K(P;n(P;v- 



The higher derivative terms in the last line can be canceled up to lower derivative terms if we chose 

1 



D-2 



ln(2K^Fit), 



(B13) 



(B14) 
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where k is an arbitrary positive constant. With this choice of the conformal factor, E^i, and E are reduced to 



E = 2«'^5^ (v^5'*''e-(^-2)-FKa,.^) + F^. (B15) 

These expressions include only up to the second derivatives. However, they include not only the conformally trans- 
formed metric g^^^ and the original scalar field (j) but also the new field oj. 

Thirdly, the new field w can be canonically normalized by defining the normalized field by 



Ki; = ^{D-1){D-2)lo = J ^-^ H2K^Fn). (B16) 



With this normalization, 



E = 2k'' 



1) ' 



(v^5''''e-^'^'^FKa..^)+F^, (B17) 



where constants a, /? and 7 are defined by (B5). Therefore, the equations of motion i^^^ = and E = ^ are given by 
(B6). 

Fourthly, it is easy to show the following equality. 

V'* (e-T«^i;^,) - le-««^i;</>,, = \E^i>.u, (B18) 

where 

E^ = V'^V^V - aKe-^-^^F - pK^e-^^^FK9^''d^(l>d,,(l) + ^e'^'^^R, (B19) 

and V is the covariant derivative compatible with the conformally transformed metric g^,y. This means that (B7) is 
derived from (B6). 

Fifthly, the variation of the action / defined by (B2) can be calculated as 



= J dPx^f^g [e-''''^Ep„gP^r"^9i.u + e'^'^'^ES^ + E^6i;] . (B20) 



Thus, the set of equations of motion derived from the action I is the same as that from I. 
Finally, by using the relations 

R = e^'^'I'R + 2-f-^K^'if^ + kV'^V' 
V>'f^ = e^'^^WV^i; - 'jKef^^^r"^,^^,., (B21) 

it is shown that i[gij.„, ip, (p] = 'l[g^l„, 4>]. 

APPENDIX C: HARMONICS ON A PLANE 

In this Appendix we give definitions of scalar, vector and tensor harmonics on an n-dimensional plane. For definitions 
and properties of more general harmonics, see Appendix B of ref. [40] and Appendix A of ref. [41] . 
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1. Scalar harmonics 

The scalar harmonics are given by 

Y = exp{-ikjx^), (CI) 

by which any function / can be expanded as 

/ = jcPk cY, (C2) 
where c is a constant depending on ki. Hereafter, we write ki as k in most cases, and sometimes we omit it. 

2. Vector harmonics 

In general, any vector field Vi can be decomposed as 

Vi = V(T)i + dif, (C3) 
where / is a function and W(T)i is a transverse vector field: 

S^^diV^T)j = 0. (C4) 
Thus, the vector field Vi can be expanded by using the scalar harmonics Y and transverse vector harmonics V^T)i 

as 

Vi = J d^k [c^T)V^T)^ + C^L)diY] . (C5) 

Here, C(t) and C(^l) are constants depending on k, and the transverse vector harmonics V(^T)i is given by 

^(T)i = Ui exp{-ikjX^), (C6) 
where the constant vector Ui satisfies the condition 

S'^kiUj = 0. (C7) 
Because of the expansion (C5), it is convenient to define longitudinal vector harmonics V(i)i by 

V^L)^ = ^^Y = -ikiY. (C8) 

3. Tensor harmonics 

In general, a symmetric second-rank tensor field Uj can be decomposed as 

ti] = t(T)ij + diVj + djVi + fSij, (C9) 
where / is a function, Vi is a vector field and t(T)ij is a transverse traceless symmetric tensor field: 

6'H^T)ij = 0, 

S"'dit^T)i'j = 0. (CIO) 

Thus, the tensor field tij can be expanded by using the scalar harmonics Y, the vector harmonics V(^T)i and V(^L)i', 
and transverse traceless tensor harmonics T(^T)ij as 

tij = j (Pk [c(^T)T(T)ij + C(^LT){diV(^T)j + djV(^T)i) 

+CiLL)idiV^L)j + djV^L)i) + C(Y)YSij] . (Cll) 
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Here, C(^t), c^lt), and C(^y) are constants depending on k, and the transverse traceless tensor harmonics T(^j')ij 

is given by 

T(T)ij = Sij exp(-iA;i/x''), (C12) 
where the constant symmetric second-rank tensor Sij satisfies the condition 

d kiSi'j = 0, 

S'hij = 0. (C13) 
Because of the expansion (Cll), it is convenient to define tensor harmonics T(^LT)ij, T(^LL)ij, and T^Y)ii 

T(LT)ij = C)iV(^T)j + C)jV(^T)i, 

= ^i(uikj + Ujki)Y, 

2 -/ , 

'^kikj ~|~ S *^ kif kjf I 
n J J J 

T(Y)ij = 5iiY. (C14) 



APPENDIX D: DETAILED ANALYSIS OF LINEARIZED GRAVITY 



In this appendix we show that linearized Einstein gravity in Minkowski background is recovered at distances longer 
than the length scale = y/cxK and at energies lower than For simplicity we consider the q — 1 case, but 

extension to a general q (> 1/2) should be straightforward. 

For this purpose we investigate in detail perturbations around the flat FRW background in the longitudinal gauge 
and take the kH — > +0 limit in the end of the calculation. This strategy makes it possible to analyze linear 
perturbations around Minkowski background, on which the action of the scalar field perturbation is not apparently 
well-defined. For simplicity, we assume that the FRW background is driven by the scalar field </>. Of course, we 
include a general matter stress-energy tensor into perturbations and carefully investigate how they couple to gravity. 
A critical point for this analysis is that constraint equations ((Oi)-componcnts and the traceless part of (zj)-componcnts 
of Einstein equation) prevent </> from fluctuating freely, and it turns out that the (j) perturbation is of sufficiently high 
order in the kH expansion that the perturbative analysis is under control. 

As stated in Sec. VII the result of the analysis is simple. The linearized gravity in our model is similar to that in 
the theory R/2k^ + aB? in the sense that the only difference is the existence of an extra scalar-type massless mode 
in our model. The extra massless mode is decoupled from the matter stress-energy tensor at the linearized level. 
Hence, as far as we are concerned with classical, linear perturbations generated by matter sources, these two theories 
give exactly the same prediction. Therefore, in our model the linearized Einstein gravity in Minkowski background is 
recovered at distances longer than the length scale /* = ^/an and at energies lower than 

Let us begin by defining the metric perturbation Sg^v by 

5m- = gf} + Sg^^., (Dl) 

where 

g^^^dx^'dx" = -dt^ + a{tfSijdx'dx^ (D2) 

{i = 1,2,3). In order to investigate hnearized gravity in Minkowski background, in the end of calculations, we shall 
take the limit 

kH +0, (D3) 
keeping d^H/H"^^ {n = 1, 2, 3) finite, where H = dta/a. In the following, we shall denote dt by an overdot. 
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1. Gauge choice and three types of perturbations 



In order to take the full advantage of the symmetry of the background spacetime, we expand the metric perturbation 
by harmonics on the 3-plane: 

Sgi^^dx^dx" = hooYdt^ + 2{h^T)oV(T)i + h(^L)oV{L)i)dtdx' 

+{h(T)T{T)ii + h(LT)T(LT)ii + h(LL)T(LL)ij + h(Y)T(Y)ij)dx'^ dx\ 

(D4) 

where Y, V(t,l) cind T(^t.lt.ll,y) ^re scalar, vector and tensor harmonics, respectively, and the coefficients hoo, 
h{T,L)Oj h(^T,LT,LL,Y) depend only on the time t. In this expression and hereafter, we omit the 3-momentum ki and 
the integration over 3-momenta. See Appendix C for definitions of the harmonics. 
We also expand the perturbation of the scalar field ^ by harmonics as 

(f) = (P^'^^t) + 6(j)Y, (D5) 

where ^^^^ (t) is the background depending only on the time t and the coefficient 5(j) of the perturbation also depends 
only on the time t. In the following we shall derive the equation of motion linearized with respect to ft,'s and S4>. 

Now, it is easily shown from the orthogonality among different kinds of harmonics that each of the following sets of 
variables form a closed set of equations: (S) /loo, ^(lo)i ^(ll), hy and S(j); (V) /i(t)Oj ^(lt); and (T) hr- Perturbations 
in different sets are completely decoupled from each other at the linearized level. The perturbations in the category 
(S) form a spin-0 representation and are called scalar perturbations. The perturbations in the category (V) form a 
spin-1 representation and are called vector perturbations. Finally, the perturbations in the category (T) form a spin-2 
representation and are called tensor perturbations. 

The metric perturbation dg^i, includes not only physical degrees of freedom but also gauge freedom. An infinitesimal 
gauge transformation is given by 

^9lJ^v ^ ^Qixv ^IJ.\v ^V\jJ,1 

5(j)Y d(f)Y - ^^'d/.cj)^'^^ , (D6) 

where is an arbitrary vector field and a semicolon denotes the covariant derivative compatible with the background 
metric gjS) ■ Hence, by expanding the vector in terms of harmonics as 

f^rfx'^ = ^oYdt + i^iT)V^T)^ + ^(L)V^L)^)dx\ (D7) 

we obtain the following infinitesimal gauge transformation for the expansion coefficients in (D4) and (D5) with ki ^ 0. 

/loo hoo — 2^0 7 
h(T)o h(T)o - o^ia^'^i.iT))' , 
h(L)o /i(L)o - 6 - a^(a"^C(L))', 

h(LT) —>■ h(^LT) - C(T)i 

2 

S(f) ^ 6(j) + ^o^^°\ (D8) 
Prom the gauge transformation (D8) it is evident that for modes with fc, ^ 0, we can choose the gauge so that 

h(L)0 = /i(LT) = = 0. (D9) 

With this condition, the gauge degrees of freedom represented by ^(t) and ^(i) are completely fixed for modes 
with ki ^ 0. On the other hand, we do not need to consider modes with fc, = since they can be absorbed into the 
background without loss of generality. 

We also take into account the stress-energy tensor T^i, of other fields, consider itself as a perturbation around 
= and expand it by harmonics as 
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Tf.^dxi'dx" = TooYdt^ + 2(T(T)oV(T)i + ^L)oV(L)i)dtdx' 

+ iT{T)T(T)tj + T{LT)T(LT)ij + T(LL)T(LL)ij + T(Y)T(Y)ij)dx'' dx^ . (DIO) 

Since gauge freedom is already completely fixed, all the coefficients r's represent gauge-inequivalent degrees of freedom. 

Hence, in the gauge (D9), all gauge-inequivalent degrees of freedom are represented by the following variables: (S) 
^00, ^'t'l ''00, '''(L)Oj '''(LL) and T(y) for scalar perturbations; (V) /i(t)Oj T(T)a and T(lt) for vector perturbations; 

and (T) h^r) and Tj^^) for tensor perturbations. As already mentioned, scalar-, vector- and tensor-type perturbations 
are completely decoupled from each other at the linearized level. Hence, in the following, we shall analyze each type 
separately. 

For vector- and tensor-type perturbations, as shown in Sec. VH, there is no difference between weak gravity in our 
model and that in Einstein theory. Hence, in the following we consider scalar-type perturbations only. 



2. Scalar perturbations 



For scalar-type metric perturbations, we have two gauge-inequivalent degrees /iqo and By introducing <i>(i) 

2; 



(Dll) 



and ^{t) by $ = — /ioo/2 and ^ = —a h(Y)/'^i the perturbed metric is written as 

g^^dx^dx" = -(1 + 2$y)dt2 (1 - 2<l!Y)a^5ijdx'dx^ , 
For this metric, the perturbed Ricci tensor is 

R^^dx^dx" = R^^^dx^'dx" + SRooYdt^ + 2dR(^L)oViL)idtdx' + {SR(^LL)T(LL)ij + 5R^Y)T(Y)ij)dx^dxK (D12) 

up to the linear order in $ and ^, where R^^v is the background Ricci tensor 



(0) 
■00 



and 



5Rqq = 2,^ + 6H^ + 3H^ 

6R^L)o = 2if + 2H^, 
1 



Rf = iH + 3H^)a^5ij, 



(D13) 



5R 



(LL) 



5R, 



{y) 



4 k' 



Ik' 



* + 6H^ + H^ + 2{H + 3H^){^ + $) + — ^ * , $ 



(D14) 



Hence, the perturbed Ricci scalar and the perturbed Einstein tensor are 
R = r(0) + 6RY, 

G^^dx'^dx'' = G^^^dx^dx'' + SGooYdt^ + 2SG^L)oV(L)idtdx' + {5G^LL)T(LL)ij + 5G(Y)T(Y)io)dx'dxi 



(D15) 



up to the linear order, where R^^^ and G^^} are the background Ricci scalar and the background Einstein tensor 



R^°^ =Q{H + 2H'^), 



'(0) 
00 

-(0) 



-{2H + 'iH'^)a^5ij, 



(D16) 



and 



5R = -2 



P 1,2 ■ 

3^ + 12iJ* -I- 3H^ + 6{H + 2if2)$ + 2—* - — $ 



<5Goo = -2 



3iJ* 



SG 



(L)0 



6G, 



(LL) 



SR 
6R 



(L)0, 



(LL), 



5G(Y) = 2a^ 



IP 



* -I- 3H^ + H^ + {2H + 3H^){^ + ^) + — ^{^ - $) 

3 ci 



(D17) 
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The Bianchi identity 

WG^^ = 0, (D18) 
where V is the covariant derivative compatible with the perturbed metric g^i,, is reduced to 

SG h'^ 

SGoo + SHSGoo + ^H—^ + ^<5G(i)o = {iG^^o + a-^S'^Gf) ^ + 2G^"o^$ - 2a-^ 5'' G^ H {-^ + $), 

<5G(,)o + ^H5G^L)o + —G^LL^ - ^ = ^ (s^S^ + a-H-G^) $ + ^a-J-cg)*. (D19) 

These two equations arc, of course, satisfied by the components shown in (D17). 
As for fields other than (j>, we have the stress energy tensor 

T^l^'^'dx^dx-' = TooYd^ + 2T(L)0ViL)idtdx' + {T^LL)T(LL)ii + ^Y)T(Y)ij)dx' dx^ . (D20) 

The conservation equation 

^ixrpother ^ P21) 

where V is the covariant derivative compatible with the perturbed metric g^j^v-, is reduced to 

Too + 3iJroo + SH — — + — T (i)o = 0, 
fiDo + 3ffT(^)o + ~nLL) - 5^ = 0. (D22) 



On the other hand, the stress energy tensor of ^ (with g = 1) is given by 
where 

L = aR^ + Lkin - V, 



(D23) 



2/ 

and a prime applied to / denotes derivative with respect to R. We now expand T^^j/ up to the linear order: 

T^f.udx'^dx" = T^^J^dx'^dx" + r^ooYdf + 2T^i^L)oV(^L)idtdx' + {T^(LL)T(LL)ij + T^{Y)T(^Y)zi)dx'dx^ . (D25) 
The background T^^J^ is given by 


V 

where 



= TT(^^^> -L- &H{LRy + 6{H' + H)Lr, 
P4>+P4> = 7r(/)(°) + 2(Lij)- - 2H{LRy + AHLr, (D27) 



i = ai?2 + ^ - y, 
LR = 2aR- ^^772-' 



2/2 

i? = 6(fl' + 2ij2). (D28) 
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and a dot represents the derivative with respect to the proper time t. The hnear part is 

1.2 

-3(H + H^)SW + W'-^^SRoo + 3H6W + —SW - 3li^(°)*, 
r4,(^L)o = F^°^^°^S(t) + W^°^SR(^L)o - SW + H6W + 



where 



2 

+2ai?W^ii + (ii + 3il2)^|^ + VF(°)o-2(5i?(y) -SW- 2HSW - -—SW 

+2VF(o) (^f + $) + 4iIW"(°) + + (^2* + , (D29) 



= 1/(0(0)), 

i^(o) = i 

~ /(i?(o))' 

(0) ^ 

(0) ^ 2[/-(it;W)]^-/(Ji;W)r(ii!W) 



= -F(^^<^W25i? + 2F%^^°^H - 2i^(°VW5.^ - AaSR. (D30) 



^S<^(°)^^i? + 2F(i 

Note that there is an identity 

4^^, W) _ 1 /o^(0) ^ „-2ry^(0)\ ^ ^ 2 ^i,-^(0),T, _ p(0) 



where i?^^'' is given by 



= - 



^(0)^(0)1 ■ _ 3ifi^(0^(0) _ 



and gives the field equation e'^'^ — for the background <p. The identity (D31) is equivalent to the (L)-coniponent of 
the linearized conservation equation of T^^^. Hence, combining this identity with the Bianchi identity (D19) and the 
conservation equation (D22) for other fields, it is concluded that the (y)-component SG(^y) = i^'^(,t^{y) + T(^y)) of the 
linearized equation of motion follows from other components. On the other hand, the 0-component of the linearized 
conservation equation of T^^^ is not an identity but is equivalent to the linearized field equation of (j). Therefore, the 
00-, {L)0- and (LL)-components give all independent equations. 

Actually, the {L)0- and (LL)-components of the equation of motion can be, in a sense, considered as constraint 
equations since we have imposed the gauge condition hf^^p = h^^LL) = 0. The (LZ/)-component SG^^ll) = i^^{t4,{ll) + 

T(LL)) is 

k'^5W = - (l - + 2kV(z,z,). (D33) 

The (L) 0-component SG(^l)q — k^{t^[l)o + '''(i)o) i^i by using the (LL)-component, reduced to 

2 

^2^(0)^(0) 5^ = (l - K^W^o)^ (^+ + + y + 3*_) + k\2^ll) - 2Ht^ll) - mo)- (D34) 
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To be precise, this equation is SG(^l)o - 2a((5G(i,L)/a)' = k'^ {t^{l)o + T(l)o) - '2'K^a[{T^(^LL) + T{LL))/a\ ■ Here, we have 
defined 



*+ = * ± $. 



(D35) 



The remaining 00-componcnt 8Gaa = t^(''",^)00 + ''00) is much more comphcated than the above two components since 
it is a dynamical equation while the above two are constraint equations. By using the {LL)- and (L)O-components, 
the 00-component is reduced to 



0, 



(D36) 



where 



Co+ = -K 



1 _ i«2^(0) _ 1^2^,(0)^(0)2 ^ SaK^ijW ] + 1^4^(0) 



^ y F(o) + ^(0) ^(0)0(0) j 

C2- = 3K2F(°)<^(o)^ 



Co- = - ai?(°)2) - in\H + i/^) (1 _ «2^(o) j ^ 3^2^,(^)^(0)2^ _ 3^4^(o) 



^ (/)(0) ^(0)0(0) j 



Soo = — K^Too + K;^f(j;,)o — 



f(ii) + 2Hf(LL) - {AH + m^)T(LL) + —r(LL) 



(D37) 



The (LL)-component (D33) can also be rewritten as a differential equation for by using the (L)O-component 
(D34) as follows. 



D2-K^^- + -Dl-K*- + Do-*- + S, 



''{LL) = 0, 



(D38) 



where 



Do+ = i/^^FW,^W2 + 



' kF^ 
F(0) ^(0) 



kH (1 - K^WW) - ^K^W^M - K [kH (1 - K^W^^^) 



2 ^(0) 



2l2 



K^k 



+ Qk'^{H + 2H'') 



Do 



3 Ac^f (") 



(^^,2^0^°^^+ 4a), 
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Dq- = 



S^LL) = -n\2f^LL) - t(L)o) + + '^C^n^D - 2Ht(ll) - T(^)o) + ^K^Hr^r^L))- - -J(^nLLy (D39) 

Now let us solve the equations (D36) and (D38) with respect to up to the zeroth order in the kH expansion. 
For this purpose, as we shall see, we actually need to keep terms up to the second order in kH in intermediate steps. 
First, let us introduce a small dimensionless parameter e, assume that 

kH = 0(e), 

^ = 0(e°), (D40) 



and estimate orders of other quantities as 



= 0(6-'*'"), 
= 0(e2), 

^2^(0)^(0) ^ o(e-i). (D41) 



Next, K^H is shown to be of order by using the equation of motion. 



'^^ = 2(1 -.^^(0)) = 

This means that kH can be considered as a constant up to the order O(e^), which is sufficient for our purpose. As a 
consequence, we can regard F^'^\ F^^ , F^j^j^ and i?^"^ as constants up to the order O(e^). The estimate (D42) leads 
to the expansion of as 

^ = ^ [1 - 2H{t - to) + 2H\t - tof + 0{e^)] , (D43) 

around the time t = to when an experiment is performed, where ao = a{to) is a constant. Here, we have assumed 
that the duration t — to of the experiment is much shorter than the cosmological time scale H~^. Of course, this 
assumption is justified since we shall take the kH +0 limit in the end of calculation. 
Thirdly, from (11) and the statement after that, we obtain the estimate that 

5!^^£!!!i!!la = 8+i = owH») = o(.>), (D44) 

C 3 6 

where we have used (D42). With this estimation and the background equation of motion, we obtain 

^ + ^ = -3.i^ ^=Oie% (D45) 

0(0) ^2^(0)0(0) ^ ' ^ ' 

By using the background equations of motion, we also obtain 
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(yW - ^ 3^2^2 (^^ ^ ^2^(0)^ ^ ^2^ ^ 2/^2^^) - i^VW - ^K^W(O) = S/^^if^ + 0(e*). (D46) 

Fourthly, let us expand all other relevant quantities as 

CO 



i=0 



Too = Tpofi], 
i=0 

oo 

i=0 
oo 

i=0 
oo 

=5ZT(ii)[i], (D47) 



i=0 



where terms with the subscript [i] is of order 0(e*). According to the expansion of r's, we obtain the following 
expansion of the conservation equation (D22): 

T'OOlO] + — '''(i)0[0] = 0, 

T(i,)o[o] + :T-2^(Li:-)[0] 12— -0 1^48) 

o Ug an 



in the order 0(e°), 



'n30[i] + ~2'^(i')o[i] + 3-ffToo[o] + 3-ff — 2 2— T(j:,)o[o]-ff(i - to) = 0, 



T{L)o[i] + o 72'^(ii)[i] - + 3-ffT(i)o[o] - ^-2nLL)[o]H{t - to) + 2-^^^H{t - to) = 0. (D49) 

O (Xn tin O tifi tin 



in the order O(e^), and 



T'oop] + — T(L)o[2] + ^Htoo[i] + 3H — 2 2— T(i)o[i]-H'(i - io) 

CI QrQ QlQ 



-6Jip.H^^t - to) + 2'-^L)o[o]H"{t - to)' = 0, 

Uq tto 

T^(I,)0[2] + ^-2T(LL)[2] 72— + '^-'^^(i)0[l] " o 72^(ii)[l]^ i* " *0j + ^— -5— -« (* " ^0 j 

O ttQ UiQ O tig ttg 

+ ^^T(^L)[o]ff'(t - to)' - 2^^H\t - to)' = 0. (D50) 

in the order O(e^). These conservation equations arc used throughout the forthcoming calculations. 
Fifthly, the coefficients of the equations (D36) and (D38) are correspondingly expanded as 

00 

00 

'S'oo = ^ 'S'oo[j] ) 
00 

Di± = y^-Dj±[»], 

i=0 

00 

S{LL) ='^S(^LL)li], (D51) 



i=0 
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where 

C2+[o] = -1, C2+[i] = 0, C2+[2] = -Aan''R^°\ 
Ci+[o] = 0) ^i+[i] = -kH, Ci+[2] = 0, 

^2^2 ^2^2 ^2^2 ^2p 

Co+[o] = 2~' ^o+[i] = 2 — ^H{t - to), C'o+p] = -2 — 2~-^^(* ~ *o)^ ~ — 2"^^^-^^°^ 



C2-[0] = ^2-[l] = C'2-[2] = C'l_[0] = Cl-[1] = C'l-[2] = C'o-p] = ^0-[l] = C'o-p] = 0, 

o 4 4* o4" o2 

"^OOpl = 7'00[0] + K '''(L)0[0] - 2k T(ii)[o] - 2 —K T(^LL){0], 

ag 

^00[1] = —1^ T00[1] + 1^ ''"(i)0[l] - -^1^ '^iLL)[l] — ^ 2~'^ 

Og 

ttg 

c 4 i4- o4" o 2 

•^OOp] = T-QOp] + ''"(i)0[2] - '''(Li)[2] - ^ 2~'^ '^(i'i')[2] 

Og 

^2^2 

+2K^iJf(i,i)o[i] - 3K''iJr(i)o[i] + 4 — 2-«^'r(iL)[i]-H'(t - ^o) - 4 — 2-K;^T(z,z,)[o]i?^(f - t^f 



a, 







«;2^(0) 

D2+[o] = -1, -D2+[i] = 0, D2+[2] = -4aK2iiW + 6a^;(Q^, 



Di+[o] = 0, -Di+[i] = -kH, -Di+[2] = 0, 

£'o+[o] = Do+[i] = 0, £>o+[2] = 2a— 2 

Da.io] = £>2-[i] = 0, = 6"— 

■Di-[o] = = -Di_[2] = 0, 



-Do-[o] 


= Do_[i] = 0, 


-Do- [2] = 


k2^^(0) /I 


S{LL)[0] 




'2k'^^ll)[o] 




%L)[1] 


= i^'^nmi] - 








= l^'^T'{L)0[2] - 


2K^f(ij:,)[2] 





K2i^(0) 



Now it is time to expand the equations (D36) and (D38). For this purpose we shall use the conservation equations 
(D48), (D49) and (D50). In the order 0{e°), by subtracting (D38) from (D36) we obtain 

^2^2 / ^2 \ 

— 2-*+[0] = '^00[0] + 2-2 T-(LL)[0] • (D53) 

"0 \ "0 / 

It is easily shown by using the conservation equations that both (D36) and (D38) are simultaneously satisfied by this 
solution in the order 0(e°). 

In the order O(e^), by subtracting (D38) from (D36), substituting the solution (D53), and using the conservation 
equations, we obtain 

^2^2 / fc^ \ 

= -K^ Too[i] + 2—T(^LL^[i] + 3iJr(z,)o[o] + 2Too[o]H{t -to)) . (D54) 

Og V flg / 

It is easily shown by using the conservation equations that both (D36) and (D38) are simultaneously satisfied by this 

solution in the order 0{e^). 

In the order O(e^), by subtracting (D38) from (D36), substituting the solutions (D53) and (D54), and using the 
conservation equations, we obtain 
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;2-*+[2] 



p(o) 



(*_[0] - 2K\LLm) - (*-[0] - 2KV(ii)[0]) 



[2], 



(D55) 



where 



"0 

-R_[2] = -K Too[2] - 2 — 2-^^ '^(ii)[2] - 2k roo[i]-H(i - io) - 3k; -H-r(i)o[i] - 6a — 2""^ 

2l2 



+ 



^2^(0) 



K'k 



i^\o[o] - Qi^\L)o[o]H^{t - to) + 8a — ^k'^R''°''t(^ll)[o]- (D56) 



By substituting this back to (D36) and using conservation equations, we obtain 



□ 



(1 - UaK^n) - 2«V(iL)[o]) + 4aK^ ( -roo[oi + 3^^) 



0. 



(D57) 



Note that ^'-[o] appeared only in the equations of order O(e^) (and higher). This is the reason why we had to keep 

terms of up to order O(e^). 

By using the conservation equation, it is easy to show that the expression roo[o] + 3r(y)[o]/ao) (^57) is 
rewritten as 



K I -T"00[0] + 3 



T(Y)[0] 



Hence, (D57) is rewritten as 



□ 



(1 - l2aK^D) *_[o] + 4aK^ ( 3*+[o] + — *+ 



-[0] 



r(LL)[0] 



0. 



(D58) 



(D59) 



Finally, by taking the kH — > limit in (D53) and (D59), we obtain the following equations for linearized gravity in 
Minkowski background. 



□ 



(D60) 



+ (too + 2fcV(ii)) = 0, 

(1 - 12aK^D) ^_ + Aan^ ^3*+ + - 2k\ll) = 0, 

where we have redefined the spatial coordinates as oqx* x* and, thus, 

□ = -df - fc2 (D61) 
The perturbed metric, after the kH +0 limit and the above redefinition of the spatial coordinates, is 

= -(1 + 2$y)rft2 + (1 _ 2^Y)6ijdx'dx^ , (D62) 
and ^± = ± The conservation equation is reduced to 



too + k r(i,)o = 0, 

4 2 

T(L)0 + 2 A T-(LL) - T(y) = 0, 



(D63) 



where r's are defined by (D20). 
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